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Q I Abstract 



The antideuteron production rate at high-energy heavy ions colhsions is calculated 

basing on the concept of d formation by antinucleons which move in the mean field 

of the fireball constituents (mainly pions). The explicit formula is presented for the 

^f) • coalescence parameter B2 in terms of deuteron binding energy and fireball volume. 
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O ' 1. Introduction. Recent measurements have reported the production of antideuterons 

r-| ■ in heavy ion colhsions P], The theoretical description of this interesting effect 

J^. is complicated, because the antideuterons are produced at intermediate stage of fireball 
OhI evolution, when the density of hadronic matter is rather low, but the particle collisions are 
still important. In other words they are produced at the "dense gas" stage of the fireball 
evolution. 

Here we present the theoretical picture of this stage and calculate the d production. The 
basic ideas of our approach are the following. The dominant mechanism of d production is 
5^ i the formation of antideuterons through the fusion reaction p + n ^ d. The fusion reaction 
is not possible if all participating particles are on mass shell. However, in the fireball at 
the "dense gas" stage of its evolution p, n, d are not on mass shell, since they interact with 
surrounding matter. The interaction with the fireball constituents leads to appearance of the 
mass shift and widths of all particles propagating in the medium (or width broadening for 
unstable ones), analogous to refraction and attenuation indices in case of photon propagation. 
The fusion reaction rate is strongly enhanced in comparison with the main process of d 
production in vacuum p + n — > d + tt. Another important ingredient of the theoretical 
picture is the balance of the deuteron formation and desintegration rates. This balance 
is achieved because of large number of produced pions and high rate of n + d collisions 
leading to d desintegration. The balance does not imply a statistical equilibrium, but rather 
a stationary process, like a balance in the isotope concentrations in a radiative chain. The 
formation rate p + n ^ d vanishes, when d size increases, i.e. its binding energy e ^ 0. This 
fact explicitly manifests itself in our calculations. The previous theoretical investigations of 
the problem were performed in statistical models [HE], in the model of (anti)nucleon sources 
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in the fireball (H] and in the Wigner function approach (see [3 |H] and references therein) . In 
all these approaches the interaction of nucleons, forming d (or d) with the fireball constituens 
as well as d (d) desintegration was not accounted (in jUini the results do not depend on e). 

According to the dominant coalescence mechanism it is convenient to characterize d 
production in heavy ion collisions by the coalescence parameter: 
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where we can put d^Np/d^pp = d^Nn/d^Pn, Pp = Pn = Pdf^- what follows, we will consider 
only the central heavy ion collisions. 

2. Theory. Consider the "dense gas" stage of fireball evolution, which follows after so- 
called "chemical freeze-out" stage [Olllllj. Assume, that particle propagations at this stage 
may be described classically using kinetic equations. We use the notation qi{x,p), i = p, n, d, 
TT ... for the double densities in coordinate and momentum spaces, nj(x) = / qi{x,p)d^p are 
the densities. {qi{x,p) - are Lorentz invariant.) Let us work in the cm. system of colliding 
ions. The kinetic equation for gj(pj, x) reads: 
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is the d 4-velocity, the ellipsis mean similar terms for collisions 



of d with other constituents of the fireball {p,n etc.) 



and v~\ v^i 

pn 



are the differences of p, 



n and vr, c? velocities Vp^ = |vp — v^l etc. The terms, when d appears in the momentum 
interval + Apj due to elastic collisions are neglected. Necessary applicability condition of 
(2) is A = p~^ -C d, where d is the mean distance between fireball constituents. 
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where Ep, En, 



Ed are p,n and d total energies, 3/4 is the spin factor and g is the coupling 
constant of low energy effective pnd intraction (in the d cm. system) . The value of was 
found by Landau jTT] from the requirement of coincidence (at the deuteron pole) of the pn 
scattering amplitude in effective theory with the amplitude in the Bethe-Peierls theory of 
the low-energy pn-scattering ^j. In the limit of zero range of nuclear forces g^ 
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1287r mNy/rrii^e, 



(4) 



where is the nucleon mass, e = 2.2 MeV is the deuteron binding energy. The account of 
non-zero range tq increases g^ by a factor of (1 — y^m/v^ro)"^ ~ 1.6 [121 US] • 

The mass of the particle moving in medium is shifted being compared with its vacuum 
value. Similarly, due to interaction with medium constutuents, the width F appears (or 
width broadening, if the particle has its proper width). The mass shift Am{E) and T{E) 
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are expressed through the forward scattering amphtude f{E) of the particle on medium 
constituent (see ^1^] and references therein). 

Am{E) = -27i — Ref{E) (5) 

T{E) = 47r-Im/(E) = ^a(E), (6) 
m m 

where E, p and m are particle energy, momentum and mass, n is the density of the con- 
stutuent in medium. Eqs. (5), (6) take place in the system, where constutuents are at rest. In 
case of moving constituents the corresponding Lorentz boost must be done. (By definition 
Am and F are Lorentz invariant, for details see [TT)j). 

Therefore, p, n and d in the reaction p + n ^ d can be considered as Breit-Wigner 
resonances with varying masses distributed according to the Breit-Wigner formula. In pro- 
cess of the fireball expansion these Breit-Wigner resonances smoothly evolve to their stable 
counterparts. So we integrate the first term in the r.h.s. of (2) after substituting (3) over 
the masses m' of the Breit-Wigner resonances: 

,,,,,, rp-/27r Tn/27T fj/27r 

dm=, dm^ drrii 



„.p, + r2/4 {mL - ninf + Tl/A (m^- - m,-)^ + r|/4 
X ^ ^ / ^ ^ q,{p,) qn{pn) S^Pp + Pn - Pi) 6{E', + E', - E',) (7) 
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where E'^ = ^p^ + m'p^ etc. We assume that the widths F <^ m are much smaller than the 
typical momenta in p, n distributions. Than the distributions qpipp) = qnipn) can be taken 
out from the integral sign at the values Pp = Pn = Pdf^- The result of calculation is given by 



Fp + F„ + F^ 2/ \ 



(the mass difference Am = mj — rrip — rrin ~ 30 MeV is small in comparison with the width 
F ~ 300 MeV and neglected in (jHJ). Later we assume Fp = F^^ = F. Fj generally is not equal 
to the antideuteron width Fj ^ 2F ~ 600 MeV. The pn system with d quantum numbers at 
high excitations will not evolve to d in the process of fireball expansion, but may decay in 
other ways. One may expect Fj < Fj. We shall keep the ratio a = Fj/Fj as free parameter 
in the calculations. However, the results weakly depend on this ratio: the variation within 
the limits < a < 1 may change the coalescence parameter ([T)) by at most times, but in 
real cases about 20%. This uncertainty is within accuracy of the whole method, estimated 
as 50%. 

The contributions of direct processes p + p d + Tc~ , n + n — >• d + 7i~^ and p + n ^ d + n'^ 
are small, alltogether about 20% in comparison with (jH}. If these processes were essential, 
the coalescence parameter B2 (HJ would be meaningless, since the antideuteron distribution 
d^Nj/d^p^ is given by complicated integral over the antinucleon distributions in this case. 

At large pp± the p spectrum decreases steeply, so the approximation pp = pj/2 becomes 
inaccurate which leads to an underestimation of B2. The case of large pp± > 1 GeV is not 
considered here. 
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Using © and performing Lorentz boost to heavy ion cm. frame, the term in square 
brackets in (j21) can be brought to the form (m^/E'j) Fj, where 



= E / d'p^ g.fe) at rest) (9) 

i 

Suppose, that the rate of antideuteron colhsions with other constituents of the fireball 
resulting in antideuteron desintegration is much larger than the rate of fireball expansion. 
This happens at collisions of heavy nuclei at high energies, when the fireball size is large 
because of large number of produced pions per nucleon. In this case one may expect the 
balance: the first term in the r.h.s. of (2) is equal to the second one and 



rj(mj/ Ej) 32 V 2 r rriTv 

The momentum distribution d^N^/d^pj entering (1) is obtained from |TU|) by integration 
over the fireball volume 



d xgj(pj,x) (11) 
Using (1), (fTUI) . pi|) and (4) (with tq correction) we find for the coalescence parameter 



Bf ='^x 1.6 Jil±^ Jd'xql{pp,x)^ ^^2^ 
^ V 2r [J d^xqp{pp,x)f 

Since the x-dependence of qpipp, x) is not known, we replace (fT^ by: 

247r3 




= =f-xl.6^/^^^^^ (13) 



where V is the fireball volume, Up and n^p are the mean and mean square p densitites in 
the fireball. (The coordinate dependence of a/T is neglected). B*^ is Lorenz invariant, as 
it should be. The volume V may be understood as a mean value of the fireball volume 
at a stage, where, on one side, hadrons are already formed, i.e., mean distances between 
them are larger than the confinement radius Rc ~ l/'^^p ~ (1/4) fm, but on the other 
side, hadron interactions are still essential. The antinucleon distributions np(r), 'n.fi(r) inside 
the fireball are nonuniform: at the dense gas stage and before it the antinucleons strongly 
annihilate in the internal part of the fireball and in much less extent in its external layer of 
the thickness of order p{n) annihilation length lann (this effect was considered in jHj). For 
this reason n^p/n^ may be remarkable larger than 1. For the same reason the antinucleons 
and antideuterons from the backside of the fireball (relative to the observer) are absorbed 
in the fireball and cannot reach the detector (see Fig. 1). Therefore, only one half of the 
fireball volume contributes to the number of registered p, n and d. The corresponding factor 
approximately equal to 2 is accounted in (fT^. 

In fact, the fireball evolution after the balance may reduce the antideuteron number 
as iVj — N^e~^^^^, where At is a typical time, required for the antideuteron to leave the 
interaction region. (Such ^ 50% reduction of K~ mesons was observed in J7j.) However, 



4 



detector 



beam 




Figure 1: Fireball at the dense gas stage. The effective volume K// is the half of the outer 
shell of thikness from which the antiprotons reach detector 



this effect does not change the coalescence parameter B2. Indeed, in this case B2 should be 
multiplied on the factor e"'^ / {e~'^^^^Y ; the time is the same both for antideuteron and 
antiproton in consideration, because they move with equal velocities. But since Fj = 2Fp 
with good accuracy (we checked it explicitly by eq ©), this factor is close to 1 regardless of 
the evolution details. 

F may be calculated if the spectrum and densities of the fireball constituents at hadronic 
gas stage are known. At high energy of heavy ion collisions (SPS, RHIC) the main con- 
ributions to Fp, F^j, Fj come from the collisions of p, n, d with pions. Therefore Fp, F^j, 
Fj essentially depend on pionic density in the fireball, the dependence of F on the densi- 
ties of other fireball constituents is much weaker. Also weak is the dependence of F on the 
spectrum of the fireball constituents, since the main contribution to Fp (Ffj, Fj) arises from 
the collisions at high energies in cm. system of p{n, d) + vr, where the cross sections are 
approximately constant. For this reason, without a serious error, for the calculation of F we 
can take the spectra from the experimental data, i.e. corresponding to the final stage of the 
fireball evolution. Moreover, since the widths enters as a/F in (|T!^ . the errors are reduced 
twice. If F is known, then by comparison with the data the parameter V~^{'n?p/n^ can be 
found, what would allow to check various models of fireball evolution. 

3. Comparison with the data. Consider the NA44 experiment at SPS (CERN): Ph + Ph 
collisions at = 17AGeV p. Antideuterons were observed at 0.6 < < 1.6 GeV and 
in the rapidity interval 1.9 to 2.1 in lab. system, which corresponds to ppt = 0.55 GeV, 
{Ep)c.m. = 1.5 GeV. The spectra and particle yields at such collisions are given in [THj . 
The number of active nucleons, participating in collision ("wounded" nucleons) iV^v and the 
number of produced pions are presented in [12]: = 362, A^^ = 1890, Q-,, = N.„/Nn = 5.2 
(see also [201 for the review of the data on heavy ion collisions). 

We accept the following model for the dense gas stage of fireball evolution [15j. (A related 
model had been suggested long ago [211 123 : it may be called Fermi-Pomeranhuk model). 
Neglect for a moment contributions of all particles except for nucleons and pions. Assume 
that any participant - nucleon or pion occupies the volume vn or v.,,, respectively. Then 

Nn n% n%Q^ 
riM = = , = — = — (14 

where n%- = 1/vn, (3 = v^^/vm- For numerical estimations we take = 0.26 fm^'^, 1.5 times 
standard nucleus density and (3 = {rT,/r^Y ^ 0.55, where r^, = 0.66 fm and rjy = 0.81 fm are 
pion and nucleon electric radii. It must be stressed, that n% is the only essential uncertain 
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parameter in our approach. Even if at dense gas stage differs from ones at the final stage, 
the arising error is essentially compensated by the appearance of Qtt both in numerator and 
denominator in (|14j) . (As was already mentioned, the nucleon contribution to F is small.) 

Check first the applicability conditions of our approach. We have: n = + rij^ ^ 
0.42 fm"^ and the mean distance between the fireball constutuents is d = = 1.3 fm. 

Evidently, the condition Xp = 1/pp <^ (i is well satisfied. The calculation of F according to 
(F = Fj/2) gives F ^ 300 MeV. (Only inelastic cross sections were accounted, the pion 
contribution comprises about 75%, the nucleon one about 25%. Note, that the value of F is 
close to the momentum integration interval in the Wigner function approach, AP ^ 200 — 
300 MeV, found in |H].) Check now the balance condition - that the probability of deuteron 
desintegration exceeds the fireball expansion rate. The former is given by 2T{miq / Ep). 
The estimation for the escape rate (or fireball expansion) is w ~ (1/4) fm~^. We have: 
2F(m7v/-Ep) ~ 2.0 fm~^ 3> 0.25 fm~^. So, this condition is also fulfilled. Even more, the 
balance condition would be fullfilled at much lower hadronic densities than chosen above, up 
to ^^5^r ~ 0.05 fm'^, i.e. up to densities not much higher, than supposed for thermal freeze-out 
[TUj . [21]. However, such low densities would lead to much lower values of i?2, than ones 
obtained in experiments. Eq.6 is legitimate, if Imf{E) -C d ^JE]. Since Im f ~ Ifm, this 
condition is not well satisfied. For this reason the value of F, presented above, has a large 
(may be 50%) uncertainty, amd, probably, is overestimated (the effect of screening). This 
fact, however, does not influence too much the value Bf^, since a/F enters ()12|) . One may 
expect, that because of their slightly larger velocities in comparison with nucleons, pions 
form a halo around the fireball. This effect also may lead to an overestimation of F. 

At the parameters used above the fireball volume comes out to be: V = 6.2 x 10^ fm^ 
(15% correction for other particles, except for pions and nucleons were accounted). This 
value is about 2 times larger, than the ones found in 19^ at chemical freeze-out and about 
2 times smaller, than at thermal freeze-out [231 121] ■ (Note, that the dense gas stage is an 
intermediate between these two.) In the case of sphere its radius is equal to R = 11.4 fm. If 
we assume, that antiprotons are mainly concentrated in the outer shell of the fireball of the 
thickness of ~ 3fm, then ti^/tt? 2 and we get for the coalescence parameter 

Bf = 3.0 X 10-^GeV^ (15) 

(We put Fj = F, or a = 1/2.) Experimentally [1], for the average value of the most central 
10% events it was found: i?2^^ = (4.4 ± 1.3) x 10""^ GeV^. However, i?2^^ strongly depends 
on centrality: the results for — 5% centrality are about 1.5 times lower. Taking in mind all 
uncertaintlies - theoretical and experimental, we believe, that the NA44 data for coalescence 
parameter are not in contradiction with theoretical expectation. 

Turn now to the STAR experiment at RHIC: Au + Au collisions at ^/s = 130 A GeV j2]. 
Antideuterons were measured at 0.5 < pt < 0.8 GeV and in the rapidity interval \Ayc.m.\ < 
0.3, 18% of central coUisioins were collected. We take Ep^^.m. = 1.05 GeV. The number of 
"wounded" nucleons in the 18% central Au+Au collisions can be estimated as Nn = 320 |2S]. 
Multiplicity of negative hadrons h (mainly, pions) was measured in [221 pseudorapidity 
7] = only and it was found an increasing of dh~ /drj |^=o by 52% comparing with the SPS 
data at i/i = 17 A GeV. But it is known that dh/drj/n=o increase faster with energy than the 
total multiplicity. We estimate = N^/Nj^ 7± 1. (A value close to the presented above, 
can be found from the data compilation [2Z1)- At Nj^ = 320 with account of 20% correction 
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for i^-mesons and hyperons V = 7.2 x 10^ fm^. The coalescence parameter is equal to 

Bl'' = 3.8 X 10-^ GeV^ (16) 

(r = 320 MeV, nV(n)2 was put to be 2). Experimentally, STAR found = (4.5 ± 0.3 ± 
1.0) X 10-^GeVl 

The main uncertainty of B2^ comes from the fireball volume V which was calculated by 
(|T^. However, the width T also depends on the fireball volume, so that B!^'^ ~ l/\/V, which 
suppresses this uncertainty twice. We expect the accuracy of our estimations (fT3j) . (fT^ to 
be about 50%. 

In E864 experiment at AGS the antideuterons were observed in Au + Pt collisions at 
y/s = 4.8y4GeV. 10% of central colhsions we selected. From the data we take: ppt = 
0.17 GeV, Ep^c.m. = 0.99 GeV. The number of "wounded" nucleons and n/N ratio are 
Nj^ = 350, Qt, = 1.6 (see ^| and references herein). In the same way as before, we 
find: = 2.8 x lO^fm^ T = 220 MeV, /^nn = 1.2 fm. In this case the validity conditions of 
our approach are at the edge of their applicability. So, the theoretical expectations for B2 
are valid only by the order of magnitude: 

5f ~ 1.5 X lO^^GeV^ (17) 

in comparison with 5f = (4.1 ± 2.9 ± 2.3) x lO^^ GeV^ 

4. Summary and Acknowledgements. The coalescence parameter B2 for the antideuteron 
production in heavy ions collisions was calculated. It was supposed, that the d production 
proceeds at the stage, when the fireball may be treated as a dense gas of interacting hadrons. 
The d production is described as the formation process p+n — > c?, where p, n, d are moving in 
the mean field of the fireball constituents (mainly pions). It was shown, that in case of large 
NT^/N]\r ratio one may expect a balance: the number of produced antideuterons is equal to 
the number of desintegrated d due to collisions with pions. The balance condition determines 
d production rate and the value of coalescence parameter B2. The later is expressed in terms 
of deuteron binding energy and mean fireball volume at this stage. The comparison with 
data demonstrates, that J-production proceeds at the stage intermediate between chemical 
and thermal freeze-out - the dense gas stage of the fireball evolution. The theoretical values 
of B2 are in satisfactory agreement with experimental data at SPS, RHIC and AGS but 
more data at various nuclei and various energies of collision and d energies would be very 
desirable. The comparison of the data with theory would allow to check various models of 
fireball evolution. 

We are thankful to G.Brown, L.McLerran, E.Shuryak for discussions and S.Kiselev, 
Yu.Kiselev, A.Smirnitsky, N.Rabin for information about experimental data. This work 
was supported in part by INTAS grant 2000-587 and RFBR grant 03-02-16209. 
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